Effect of A-site size difference on polar behavior in MBiScNbC>6, (M=Na, K and Rb): 

Density functional calculations 
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We investigate the effect of A-site size differences in the double perovskites BiScC>3-MNb03 
(M=Na, K and Rb) using first-principles calculations. We find that the polarization of these ma- 
terials is 70~90 /iC/cm 2 along the rhombohedral direction. The main contribution to the high 
polarization comes from large off-centerings of Bi ions, which are strongly enhanced by the sup- 
pression of octahedral tilts as the M ion size increases. A high Born effective charge of Nb also 
contributes to the polarization and this contribution is also enhanced by increasing the M ion size. 

PACS numbers: 77.84.Ek,71.20.Ps 



INTRODUCTION 



Ferroelectric perovskite oxides have attracted atten- 
tion for many years because of their many applications 
in technology, including memory devices, electronic com- 
ponents, piezoelectric sensors and actuators, and pyro- 
electric devicespi Besides the technological importance, 
there is also a great scientific interest in understanding 
the underlying physics. This is stimulated by recent dis- 
coveries of entirely new ferroelectrics particularly in thin 
film form^— The key to these materials is competition 
between different interactions, leading to structural in- 
stability. Ferroelectricity in perovskite ABO% oxides de- 
pends on the delicate balance of the long-range Coulomb 
energy, which favors a ferroelectric distortion, and the 
short-range ionic repulsion, which stabilizes the cubic 
phasei££ More generally, instabilities in perovskites are 
often discussed in terms of competition between preferred 
A-0 and B-0 bond lengths, which cannot both be accom- 
modated in the ideal cubic structure for arbitrary A and 
B ions. This is characterized by a so-called tolerance 
factor, t, that describes the degree of non-satisfaction 
of ideal bond lengths in the cubic perovskite structure. 
The balance between the long-range Coulomb energy and 
short-range ionic repulsion can be modified by smaller 
effects related to the A-0 hybridization, in particular 
the so-called lone pair physics associated with stereo- 
chemically active ions, particularly Pb 2+ and Bi 3+ , as 
in the technologically important Pb(Zr,Ti)03 system and 
BiFeOa . In thin films and superlattices the constraint im- 
posed by lattice matching imposes strains that push bond 
lengths away from their ideal values, leading to highly en- 
hanced ferroelectric instabilities in many cases^i!^ 

The piezoelectric performance in perovskite ferro- 
electrics is associated with polarization rotation near 
morphotropic phase boundaries (MPB). 8 ~ 13 That is, in 
materials where two different ferroelectric states with dif- 
ferent lattice strains are close in energy, an electric field 
may rotate the polarization and therefore the ferroelec- 
tric state from the zero field equilibrium configuration 
towards the other nearby configuration, with resulting 



coupling to the lattice strain. Strong ferroelectricity is 
highly desirable for piezoelectricity, but epitaxial con- 
straints on the lattice strain as in most thin films are 
not. It is of interest to examine strategies for producing 
structural frustration and perhaps enhance ferroelectric- 
ity in bulk materials. 

In this regard, recent first principles calculations sug- 
gest that A-site off-centering can be enhanced by the use 
of ions with different ionic radii at the A-site There is 
experimental support for this in the (Ba,Ca)Ti0 3 solid 
solution, which shows a strongly stabilized ferroelectric 
tetragonal phase over a wide composition range as com- 
pared to what would be expected from a classical averag- 
ing of the end-point properties. 15 However, while this is 
shown experimentally to be effective in enhancing ferro- 
electricity in (Ba,Ca)TiC>3 it is not effective in the closely 
related (Ba,Ca)ZrC>3 system.— We recently reported re- 
sults of a first principles study of the double perovskites, 
BiPbZnNbC>6 and BiSrZnNbOe, which have a large size 
difference on the A-siteJ^L We found high polarization 
values of ^80 /iC/cm 2 in these compounds, which is com- 
parable to BiFe03 (90-100 ^iC/cm 2 )^ even though these 
double perovskites have only half the A-sites occupied 
by the stereochemically active Bi 3+ ion. The high po- 
larization in these materials results comes from large off- 
centerings of Bi ions, associated in part with the A-site 
size mismatch, cooperating with moderate off-centerings 
of Nb ions that have Born effective charges highly en- 
hanced by covalency. 

Here we use density functional calculations to study 
the double perovskites MBiScNb0 6 (M = Na, K and 
Rb) in order to further investigate the effect of A-site 
size differences on the polar behavior of perovskites. In 
relation to BiSrZnNbOg, the compounds that we study 
here have lower ionic charges at the A-site and higher 
ionic charges at the B-site. The M ions have increas- 
ing ionic radii (r Na + < rx+ < ?Rb+)> but are otherwise 
chemically similar, which is helpful in understanding how 
increasing A-site size difference affects the ferroelectric- 
ity. Finally, we note that these materials could poten- 
tially be of practical use. It is known that the end-point 
compound BiSc03 (BS) strongly favors a rhombohedral 
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state over a tetragonal state^^— as does BiSrZnNbOg. 
As such these systems, when alloyed with a tetragonal 
ferroelectric can show MPBs and possibly form useful 
piczoelectrics J£ 



to chemical ordering. Thus, any off-centering is a conse- 
quence of lattice instabilities and not an artifact due to 
cation ordering in the supercell. 



II. APPROACH 

We investigated the lattice instabilities of three per- 
ovskite oxides (Na,Bi)(Sc,Nb)0 6 , (K,Bi)(Sc,Nb)0 6 and 
(Rb,Bi)(Sc,Nb)0 6 using density functional theory (DFT) 
within the local-density approximation (LDA). All cal- 
culations, including the structure relaxations, polariza- 
tions and Born effective charges, were done using a 
plane wave basis as implemented in Quantum Espresso 
packaged We used ultrasoft pseudopotentials 23 and con- 
verged our results with respect to the cut-off energy and 
Brillouin zone sampling, similar to our previous study for 
BiPbZnNbOg and BiSrZnNb0 6 .— The polarization was 
calculated using the Berry's phase method. 

There is a large charge difference between Sc 3+ and 
Nb 5+ , and this is expected lead to a strong ordering ten- 
dency on the £?-site. This is actually the case in the re- 
laxor compound PbSc^Nb^Oa (PSN), where Sc and 
Nb have rock-salt ordering] 24 ' 25 Here we consider ma- 
terial with this double perovskite type B-site ordering. 
In addition, while ordering at the A-site is less clear, 
there is a large charge difference at the A-site, which 
may lead to an ordering tendency at least short range, as 
well. There have been reports of short-range cation or- 
dering of Na+ or K+ and Bi 3+ in Na 1/2 Bi 1/2 Ti0 3 (NBT) 
and K 1/2 Bi 1/2 Ti0 3 (KBT) which support this view. 
Here we consider the ordered system, both because it is 
simpler, which facilitates the calculations and analysis, 
and because it helps to avoid computational artifacts, 
such as polarization because of the particular symmetry 
breaking chemical ordering in a given supercell. 

There are various factors at play that may drive 
ferroelectricity in NaBiScNbOg (NBSN), KBiScNb0 6 
(KBSN) and RbBiScNbOg (RBSN). These materials 
have Nb 5+ at the -B-site which favors ferroelectricity, as 
in KNbC>3. The end-point perovskite BiScC>3 has a very 
small tolerance factor t = 0.85 and a stereochemically ac- 
tive Bi 3+ at the A-site, two factors that are considered to 
be important in A-site driven materials, leading to strong 
lattice instabilities. Furthermore, the large size differ- 
ence at A-site may also enhance the off-centering of the 
A-site ions through frustration of tilt modes. However, 
in addition to A-site off-centering, low tolerance factors 
also strongly favor rotation of the -BOg octahedra. These 
rotations compete with ferroelectricity in general. In or- 
der to sort this out it is important to use supercells that 
have even numbers of units along the [100], [010] and 
[001] directions to accommodate the various Glazer tilt 
patterns^ We used 40 atom 2x2x2 supercells and as- 
sumed both A- and £?-site rocksalt chemical ordering. In 
addition to allowing various Glazer tilt patterns, these 
supercells are cubic and therefore non-polar as regards 



III. RESULTS AND DISCUSSION 

To our knowledge, syntheses of NaBiScNbOg, 
KBiScNbOg and RbBiScNbOg have not been reported. 
We therefore used theoretical LDA lattice parameters. 
We note that in oxides with heavy elements the LDA 
typically leads to a small underestimation of the lattice 
parameter by ~ 1-2%. If this is the case here, it will 
work against ferroelectricity and in favor of competing 
tilt instabilities^ We assume a pseudocubic structure to 
evaluate the lattice parameters and relaxed all atomic 
coordinates with no symmetry constraints. This was 
done as a function of the lattice parameter to find the 
equilibrium volumes for these materials. The calculated 
LDA lattice parameters for the minimum energy were 
3.98 A, 4.03 A and 4.07 A for NaBiScNbOg, KBiScNbOg 
and RbBiScNbOg , respectively. These show the expected 
trend with increase in the A-site ion size. 
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(c) RbBiScNbOe 

FIG. 1: (color online) Projections onto a (010) plane of 
the structures of the relaxed 40 atom pseudocubic (a) 
NaBiScNbOe, (b) KBiScNbOe and (c) RbBiScNbOe super- 
cells. The B-0 bonds are shown. The Bi ions are denoted by 
light purple spheres, the Sc ions by blue, the Nb ions by light 
green, and the O ions by red. The Nb, K and Rb ions are 
shown by dark yellow spheres in each panel. The O octahedra 
for Sc are shown by (semi-transparent) blue boxes and those 
for Nb by (semi-transparent) green. The solid lines show the 
unit cell boundaries of each material. 



Fig.Q]shows the structures of the relaxed 40 atom pseu- 
docubic NaBiScNbOg, KBiScNbOg and RbBiScNbOg su- 
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percells at the equilibrium volume. As expected, the 
octahedral tilting is reduced with larger A-site ion size 
difference (note that the ionic radii of A-site ions are 
r Bi 3+ = 1.31 A, r Na + = 1.53 A, r K + = 1.78 A and 
r Rb + = 1.86 A, respectively). In particular, there is rela- 
tively little tilt in RbBiScNbOe, which has the largest A- 
site size difference. Fig. [T] also shows large displacements 
of the Bi ions from the ideal positions at the centers of 
the O cages. 




Bil Bi2 Bi3 Bi4NalNa2Na3Na4Scl Sc2 Sc3 Sc4NblNb2Nb3Nb4 

FIG. 2: (color online) Cation off-centerings along the Carte- 
sian directions and their magnitude with respect to the center 
of their O cages (the 12 nearest O ions for the A-site and the 
6 nearest for the B-site) in the relaxed 40 atom NaBiScNbOe 
supercell. 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




BNSN 
BKSN 
BRSN 



LiM 



FIG. 3: (color online) Average magnitude of cation off- 
centerings with respect to the center of their O cages in the 
relaxed 40 atom NaBiScNbOe, KBiScNbOe and RbBiScNbOe 
supercells. 

The cation off-centerings in NaBiScNbOe from the cen- 
ter of their O cages is shown in Fig. [3] All cations off- 
center along [111] direction and the magnitude of dis- 
placement is almost the same for a given cation. The 
Bi ions have a large off-centering, which is a conse- 
quence of the difference in size compared to the other 
A-s\te ion and also its stereochemical activity. In addi- 
tion to the off-centering of A-site cations, there are dis- 
placements of B-site cations as well. For NaBiScNbOe, 
the average magnitudes of the cation off-centerings are 



0.62 A for Bi 3 +, 0.28 A for Na+, 0.18 A for Sc 3 + and 
0.24 A for Nb 5+ , respectively. This cooperation of all 
cations in the ferroelectric distortion is not unusual and 
is also seen in KNb0 3 , BaTi0 3 , PbTi0 3 , BiPbZnNb0 6 
and BiSrZnNb0 6 >iii^i^ The cation off-centerings in 
KBiScNbOe and RbBiScNbOe are also collinear along 
[111] direction. 

As mentioned, the LDA lattice parameters of these su- 
percells may be underestimated by 1-2%. To study the 
robustness of ferroelectric displacements, we performed 
relaxation calculations for larger values of lattice param- 
eters that lie between the equilibrium lattice parameters 
(3.98 A, for NaBiScNbOg, 4.03 A, for KBiScNbOe and 
4.07 A, for NaBiScNbOe) and the higher value of 4.23 
A. We found ferroelectric cation off-centerings for each 
such calculation, confirming the robustness of ferroelec- 
tric distortions against LDA volume errors. 

Fig. [3] shows the average magnitudes of the 
cation off-centerings in NaBiScNbOe, KBiScNbOe and 
RbBiScNbOe . The displacements of Bi and Nb increase 
as the A-site size difference is increased. The reason 
for this increase in displacement is two-fold. First, the 
larger size difference decreases the tilting of BOq octa- 
hedra and provides more volume for displacement of Bi 
and Nb ions. Secondly, the lattice parameter increases as 
the size of M ion increases. The ferroelectricity becomes 
stronger with increasing volume, as is usual in perovskite 
oxides. 34 Competing with this is the fact that, since po- 
larization is dipole per unit volume, larger unit cell vol- 
ume works against high polarization. It is interesting to 
note that in contrast to Bi and Nb ions, the off-centerings 
become smaller for the M ions, while there is no appre- 
ciable change in the displacement of Sc ion. As discussed 
below, there is little covalency between the unoccupied 
states of the M ions with the occupied O 2p states in 
these materials. This works against off-centering. Fur- 
thermore, r Na + < t*k+ < ?*Rb+ > an d as the radius of 
an atom increases, there is less volume available for off- 
centering. This explains the smaller displacement of Rb 
ions as compared to Na ions. 

We now consider the electronic structure of these ma- 
terials as it pertains to ferroelectricity. Fig. @] shows 
the electronic density of states (DOS) of the pseudocu- 
bic relaxed NaBiScNbOe and the projections onto Bi p, 
Sc d, Nb d and O p pseudo-wavefunctions (the DOS of 
KBiScNbOe and RbBiScNbOe are almost identical and 
are not shown here). The O 2p states provide the DOS 
from -5 eV to the valence band maximum and there is a 
gap that separates the O 2p derived states from the un- 
occupied metal states. The calculated LDA band gap is 
3.2 eV, 3.1 eV and 3.1 eV for NaBiScNbOe, KBiScNbOe 
and RbBiScNbOe, respectively. The actual band gaps 
are likely to be larger than these calculated ones due to 
LDA band gap underestimation. In any case, based on 
the calculated gaps, these materials are expected to be 
good insulators at ambient temperature. The O 2p de- 
rived valence bands contain significant admixtures of the 
nominally unoccupied metal states, which are mainly of 
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FIG. 4: (color online) Electronic DOS of the pseudocubic 
relaxed NaBiScNbC>6 and the Bi p, Sc d, Nb d and O p pro- 
jections. The valence band maximum is taken as the energy 
zero. 



TABLE I: Comparison of properties of NaBiScNbOe, 
KBiScNbOe and RbBiScNbOe, where r M + denotes the ion 
radii of M ion, a denotes the pseudocubic lattice parameter, 
z* denotes the average Born effective charge, S denotes the 
average cation off-centering with respect to the O cages, and 
P is the magnitude of polarization. 





NaBiScNbOe 


KBiScNbOs 


RbBiScNbOs 


r M+ (Na/K/Rb) 


1.53 A 


1.78 A 


1.86 A 


a(LDA) 


3.98 A 


4.03 A 


4.07 A 


*(Bi) 


0.62 A 


0.67 A 


0.77 A 


**(Bi) 


4.38 


4.21 


3.73 


5(Na/K/Rb) 


0.28 A 


0.22 A 


0.20 A 


z*(Na.K/Rb) 


1.14 


1.17 


1.22 


<5(Sc) 


0.18 A 


0.16 A 


0.16 A 


z*(Sc) 


4.73 


4.81 


4.92 


<5(Nb) 


0.24 A 


0.27 A 


0.31 A 


z*(Nb) 


6.42 


6.23 


5.80 


P(lll) 


73 fiC/cm 2 


77 nC/cm 2 


87 nC/cm 2 



Bi 6p, Sc 3d and Nb Ad character. This mixing of oc- 
cupied O 2p states with unoccupied d states of metals 
implies chemical bonding and is responsible for enhanced 
Born effective charges and ferroelectric off-centering, as 
in PbTi0 3 and BaTi0 3 ^ Nb 4d states are dominant 
near the conduction band minimum and this low energy 
favors hybridization with O 2p states, as discussed in the 
case of PbMg!/ 3 Nb2/ 3 3 and KNb0 3 in comparison with 
KTa03i 36 i 37 However, the Sc 3d states lie further above 
the conduction band minimum and this results in a some- 
what weaker hybridization with O 2p states. We do not 
find significant mixing of unoccupied Na/K/Rb s states 
with the occupied O 2p states as may be expected consid- 
ering the very large electronegativity difference between 
these elements and O. In NBSN, the influence of cova- 
lency is reflected in the enhanced Born effective charges 
of 4.38 for Bi, 4.73 for Sc and 6.42 for Nb as compared 
to the very modest enhancement of (z*=1.14) for Na. 
Similar numbers are obtained for KBSN and RBSN (see 
Table U. 

In contrast to the trend in the cation off-centerings, the 
Born effective charges of Bi and Nb decrease, while those 
of M and Sc ions increase when the A-site size difference 
becomes larger. These decreases in Born charges of Bi 
and Nb ions would mean lower values of polarization if 
the Born charges were independent of atom position and 
the off-centering were of the same magnitude in all com- 
pounds. In fact, however, the polarizations as calculated 
with the Berry phase method are 73 /iC/cm 2 for NBSN, 
77 nC/cm 2 for KBSN, and 87 ^C/cm 2 . 

As mentioned, high piezoelectric performance is in- 
duced by the polarization rotation from the rhombohc- 
dral state to tetragonal state near the MPB. We studied 
the tendency towards tetragonality by calculating the to- 
tal energies of NBSN, KBSN, and RBSN as a function 
of imposed tetragonal strain (0.94 < c/a < 1.06) while 
keeping the volumes equal to their equilibrium volume 
for the pseudocubic cells. In all materials, the lowest en- 



ergy is at c/a — 1.00 and the cation off-centerings are 
collinear and along [111] direction. The cations do not 
off center along [001] directions even with large c/a and 
the off-centerings remain practically along the [If 1] direc- 
tion over the range of c/a values that we tried, implying 
robustness of the rhombohedral ferroelectricity. There- 
fore, these are rhombohedral ferroelectric materials, and 
as such it is of interest to alloy with tetragonal ferro- 
electrics in order to produce a potential MPB. 

Related to this, a similar compound Nao.sBio.5Ti0 3 
(NBT) is also a rhombohedral ferroelectrio 3 ^ with Tq — 
593.15 K and polarization 38 /iC/cm 2 ^ Solid solutions 
of NBT with tetragonal BaTi0 3 (NBBT) show a useful 
MPB with polarization 20 /iC/cm 2 ^ 9 While the mate- 
rial has not proved superior to other compounds such 
as Pb(Zr,Ti)0 3 in transducers, it does have a lower free 
permittivity and a high electromechanical coupling fac- 
tor that make it an interesting lead-free actuator mate- 
rial. The materials in the present study are similar in 
that they are rhombohedral ferroelectrics and may be 
considered to be related to NBT with Ti replaced by 
Sco.5Nbo.5- By analogy, solid solutions of these materi- 
als with BaTi0 3 may provide materials that are useful 
as lead-free piezoelectrics. Another possibility may be al- 
loying with the super tetragonal material Bi(Zn,Ti)0 3 .— 

IV. CONCLUSIONS 

We investigated the polar behavior of double per- 
ovskites NaBiScNbOg, KBiScNbOg and RbBiScNb0 6 
using density functional calculations on small supercells. 
We find high polarization along [111] directions. The 
main contribution to the polarization comes from large 
displacements of Bi and more modest displacements of 
Nb ions, which, however, have high Born charges. We 
find that the polarization increases as the A-site size dif- 
ference becomes larger. These materials also have sub- 
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stantial band gaps and the rhombohedral ferroelectricity 
is robust. It will be of interest to investigate whether 
solid solutions of these materials with tetragonal ferro- 
electrics show MPBs possibly providing another avenue 
for development of lead-free piezoelectrics. 
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